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PREFACE 


This  report  discusses  the  effects  of  gravity  and  of  the 
atmosphere's  nonuniform  density  on  sound  propagation.  It  is 
part  of  RAND'  s  continuing  study  of  the  atmospheric  waves 
generated  by  nuclear  explosions. 

This  report  should  be  of  interest  to  agencies  and  contractors 
concerned  with  the  detection  of  nuclear  explosions.  It  should 
also  be  of  interest  to  other  research  workers  who  are  studying 
low  frequency  sound  propagation  in  the  atmosphere. 


SUMMARY 


A  theoretical  discussion  is  presented  of  the  influence  of  gravity 
on  sound  propagation  from  a  small  source  in  an  isothermal  atmosphere 
where  ambient  pressure  and  density  decrease  exponentially  with  height. 

A  solution  for  the  free-space  case  is  derived  which  indicates  that 
waves  with  angular  frequency  w  between  (7- l)^2g/c  and  (r/2)g/c  will 
not  be  propagated,  while  those  with  u)  between  0  and  (7-l)^"^(g/c)  cos  0 
will  not  be  propagated  in  a  direction  making  an  angle  of  9  with 
the  vertical  axis.  A  formal  solution  incorporating  appropriate 
boundary  conditions  at  the  ground  is  derived  and  discussed.  The 
field  along  the  vertical  line  passing  through  the  source  is  found 
explicitly.  A  consideration  of  the  energy  intensity  shows  that  no 
energy  is  propagated  within  a  cone  above  and  below  the  source  if 
cj  <  (7~l)^^g/c.  A  calculation  of  the  intensity  for  the  case  when 
(7-l)^2g/c  <  u)  <  ( 7/2  )g/ c  indicates  that  the  energy  flowing  from  the 
source  tends  to  concentrate  in  the  lowest  layers  of  the  atmosphere. 

The  field  for  large  horizontal  distances  appears  as  a  sum  of  a  direct 
wave,  a  reflected  wave,  and  a  surface  wave.  Reflection  coefficients 
axe  derived  and  the  criteria  for  the  surface  wave  to  be  dominant  are 
discussed. 
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I.  INTRODUCTION 


The  effects  of  the  earth1 s  gravitational  force  on  the  propagation 
of  infrasonic  waves  to  great  distances  has  been  extensively  studied  by 
a  number  of  workers/1  Interest  in  this  subject  dates  back  to  1883; 

when  the  eruption  of  the  Krakatoa  volcano  generated  a  pressure  pulse 
which  was  detected  in  widely  spaced  regions  of  the  world.  The  detection 
of  waves  from  the  great  Siberian  meteorite  in  1908  and,  in  recent  years, 

from  the  detonation  of  hydrogen  bombs  created  additional  interest. 

( 1-5) 

Early  theoretical  studies'  "  showed  that  one  of  the  principal 
effects  of  the  earth*  s  gravitational  attraction  was  that  it  permitted 
the  existence  of  horizontally  propagating  modes  trapped  in  the  lower 
atmosphere.  These  modes  were  found  to  exist  when  the  assumed  temperature- 
height  profile  of  the  atmosphere  had  no  temperature  minimum,  and  when 
the  conventional  acoustic  theories  (i.e.,  those  neglecting  the  influence 
of  gravity)  precluded  the  existence  of  such  modes. 

In  this  paper  we  shall  study  the  manner  in  which  these  trapped 
modes  (or  gravity  modes)  evolve  from  a  small  source  located  above  the 
ground.  In  particular,  we  shall  derive  expressions  for  the  acoustic 
field  in  regions  relatively  close  to  the  source  and  shall  study  the 
effects  of  the  ground  on  this  field. 

The  model  we  shall  study  is  that  of  an  isothermal  atmosphere 

bounded  by  a  flat  earth.  Although  this  model  is  an  oversimplification 

of  the  real  atmosphere  and,  as  is  well  known,  it  does  not  lead  to  a 

correct  prediction  of  the  form  of  the  pressure  pulse  which  would  be 

observed  at  large  distances,  we  feel  that  a  study  of  this  model  may 

lead  to  a  better  understanding  of  the  origin  of  the  phenomena  predicted 

at  great  distances  by  more  sophisticated  models.  In  particular,  the 

theory  developed  here  should  complement  the  recent  asymptotic  theories 

(7)  ( 8) 

of  Hunt,  Palmer,  and  Penney'  '  and  of  Weston. 

The  subject  of  wave  propagation  in  an  isothermal  atmosphere  has 
been  considered  previously  by  Lamb/1)  by  Pekeris/^  and  by  Sretenskii/11) 
Lftrajh  discussed  the  problem  of  the  vertical  and  horizontal  propagation 
of  plane  waves  and  Pekeris  extended  the  discussion  of  horizontal 
propagation  to  illustrate  his  theory  of  the  excitation  of  gravity 
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modes  by  a  small  source  on  the  ground  in  an  atmosphere  having  a 
constant  lapse  rate  in  the  troposphere.  Sretenskii  considered  the 
propagation  of  waves  from  a  point  source  in  an  isothermal  atmosphere 
with  the  influence  of  the  earth's  surface  neglected  (i.e.,  the  source 
was  sufficiently  far  above  the  earth  that  reflections  could  be  dis¬ 
regarded).  It  appears,  however,  that  Sretenskii' s  theory  is  in  error. 
His  initial  equation, 


for  the  velocity  potential  is  inconsistent  with  the  equations  of 
hydrodynamics  for  air.  Lamk/^  has  shown  that  a  velocity  potential 
does  not  exist  for  the  Isothermal  atmosphere  unless  7*1*  If 
Sretenskii's  theory  were  to  have  any  application  to  the  real  atmosphere, 
he  would  have  to  take  y  *  1.4.  Furthermore,  Sretenskii  imposed  the 
requirement  that 


3* 

3z 


0 


(1.2) 


everywhere  in  the  horizontal  plane  of  the  source  except  at  the 
source.  The  only  justification  given  by  Sretenskii  for  this  require¬ 
ment  is  that  it  is  applicable  when  gravity  is  neglected.  The  early 
portion  (Sec.  Ill)  of  this  paper  will  be  concerned  with  the  presenta¬ 
tion  of  the  correct  solution  of  this  problem. 
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II  .  MATHEMATICAL  FORMULATION  OF  THE  PROBLEM 

For  an  isothermal  atmosphere,  the  characteristic  length  is  the 
scale  height, 

H  *  c2/(7g)  (2.1) 

where  c  is  the  (constant)  speed  of  sound,  g  is  the  acceleration  of 

gravity,  and  y  is  the  specific  heat  ratio,  which  for  air  may  be  taken 

as  1.40.  The  characteristic  time  is  h/c.  We  shall  accordingly 

develop  our  theory  in  a  system  of  units  in  which  distance  is  in  units 

of  H  and  time  is  in  units  of  h/c.  (Typical  values  of  H  and  h/c  for 

the  real  atmosphere  are  8  km  and  24  sec.) 

Let  the  pressure  and  density  at  a  pointy  at  time  t  be  represented 

by  p  +■  p  and  p  +  p,  where  p  and  p  are  their  ambient  values.  For 
*  o  o  o  o 

an  isothermal  atmosphere,  pQ  and  pQ  are  given  by  the  expressions 

Po  =  c2po/r  =  Pooe-z  ,  (2.2) 

where  pqq  is  the  pressure  at  the  ground  and  z  denotes  the  height 
above  the  ground  (in  units  of  H).  The  quantities  p,  p,  and  the 
particle  velocity  are  presumed  to  obey  the  linearized  equations 
of  hydrodynamics.  For  an  isothermal  atmosphere,  with  distances  in 
units  of  H  and  times  in  units  of  H/c,  these  have  the  form: 

+  ;$>  +^zc2p  h 

cdp/dt  +  J  *  (  PQv) 

2 

dp/dt  -  c  dp/dt 


*  ~  (ecpgX)  (2.5a) 

«  0  (2«3t>) 

+  e7k\‘  (P  v)  -  0  (2.3c) 
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where  denotes  the  unit  vector  in  the  z-direction  and 

A2  -  (7  -  1 )/y2  M  .204  .  (2.4) 

These  equations  may  he  derived  in  the  manner  outlined  hy  Pridmore- 

(12) 

Brown.  The  term  in  braces  in  (2.3a)  is  included  in  order  that 

we  may  qualitatively  take  frictional  forces  into  account.  The 
parameter  e  is  a  positive  constant  which  we  shall  assume  to  be  very 
small.  The  inclusion  of  this  term  is  a  well  known  device^^  which 
is  generally  attributed  to  Rayleigh. 

The  Fourier  transforms  p,  p,  and  v  of  p,  p,  and  v,  where 

S/A 

p(u>)  =  ( 2jt ) _1  Jp(t)  eilJt  dt  ,  (2.5) 

0 

etc.,  will  satisfy  three  simultaneous  differential  equations  which 
may  be  formally  obtained  by  replacing  3/&t  by  -ito  in  equations  (2.3). 
(Here  u>  represents  the  angular  frequency  in  units  of  c/H.)  From 
these  equations  one  may  derive  a  single  differential  equation  for 
the  Fourier  transform  of  any  one  of  a  number  of  quantities  associated 
with  the  acoustic  field.  (Lamb,^  for  example,  derived  a  differential 
equation  for  the  vorticity.)  We  choose  to  concentrate  on  the  pressure. 
If  we  define 


a  z/2 
pe  ' 


(2.6) 


then  this  quantity  P  will  satisfy  the  following  differential  equation 


7;(7P)  +  (P2/u2)P  =  0 


(2.7) 


where 


e2  -  n2  -  i  , 


2  2 

cr  -  a  , 


(2.8) 

(2.9) 
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and 


if 


u>(o)  +  ic) 


(2.10) 


The  operator  7  has  components  (l/fi  )d/dx,  (l/Q  ;d/dy,  (l/pt  )d/dz  in 
Cartesian  coordinates*  Thus 


2  2  2 

~  \  2  +  ^  2'  +  2  2* 
0  dx  By  p, 


(2.11) 


The  quantity  P  plays  the  role  of  a  potential.  Once  it  is  known, 

/N  A 

p  and^v  may  be  obtained  by  use  of  the  following  formulas: 

C2p  =  e'z/2  {o2  -  (7/2 )A2  +  7A20/az)}  p/n2  ,  (2.12) 

i(u  +  l«)cP(£  =  e"z/2|Q2£p  +  (02/u,2^ezB2p|  ,  (2.13) 


where 


B  =  (1/7)  -  (1/2) 


(2.1k) 


Since  the  vertical  component  of  the  particle  velocity  must 
vanish  at  the  ground,  P  must  satisfy  the  boundary  condition 

dp/az  +  b2p  =  0  (2.15) 

at  z*0.  This  follows  directly  from  (2.13). 

In  addition,  P  must  conform  to  causality.  In  conventional 
acoustic  theories,  this  would  be  interpreted  as  requiring  that  P 
behave  as  an  outgoing  wave  in  regions  remote  from  the  source 
(radiation  condition) .  A  general  formulation  of  this  requirement 
is  that,  for  all  such  that  z  >  0,  the  quantity  Pe~iurt°,  where  t 
corresponds  to  any  time  before  the  source  is  initially  excited, 
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should  be  analytic  everywhere  in  the  upper  half  of  the  complex 
to-plane  and  should  vanish  as  Im(w)  tends  to  infinity.  This  implies, 
in  particular,  that  P  should  have  no  poles  or  branch  lines  in  the 
region  for  which  Im(a>)  >  0.  That  this  requirement  guarantees 
causality  is  clear,  since  the  inverse  transform 

p(r>  t)  =  e_Z/2  J  P(r,w)e_iwtdw  (2.16) 

will  then  vanish  for  all  t  <  t 
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III,  THE  FREE-SPACE  SOLUTION  FOR  A  SMALL  SOURCE 


To  examine  the  nature  of  the  field  close  to  the  source,  we  shall 
solve  the  free-space  problem.  The  presence  of  the  ground  is  neglected 
and  the  boundary  condition  (2.15)  is  discarded.  The  causality  condition 
is  required  for  all  rather  than  for  Just  those  for  which  2  >  0. 

A  Green*  s  function  G(^-jfo,w)  may  readily  be  found  which  is  finite 
everywhere  except  at  a  point  r  (assuming  e  >  0),  and  which  conforms 
to  causality  in  the  sense  described  above.  We  may  choose  the  normaliza¬ 
tion  of  the  Green* s  function  such  that  it  satisfies  the  equation. 

v-7G  +  O2/u2)0  =  -4n8(r-i^)/(Q2u).  (5-1) 


It  is  clear  that  this  equation  is  formally  equivalent  to  that  for  the 
ikR  / 

Green*  s  function  e  /R  for  the  scalar  Helmhltz  equation.  One  need 

only  replace  Ox,  Qy,  and  \).z  by  x* ,  y* ,  and  z*  to  cast  it  in  a  form 

(14) 

identical  to  that  given  by  Morse  and  Feshbach.  7  Thus  we  may  write 


G(r-r  ,to) 
^o 


R* 


-1  ei(p/u)R* 


(3-2) 


where 


R'  =  jn2(x-xQ)2  +  n2(y-yo)2  +  u2(z-zo)2J  .  (30^ 

The  requirement  of  causality  may  be  satisfied  by  a  suitable 
definition  of  the  phases  of  p,  n,  and  R'  in  the  complex  w-plane. 

The  branch  lines  for  these  quantities  are  taken  as  extending  vertically 
downwards  from  their  respective  branch  points--all  of  which  lie  slightly 
below  the  real  axis.  One  requires  that  the  phases  of  p,  v,  and  R'  be 
continuous  everywhere  except  at  their  branch  lines  and  that ,  for 
Im(to)  >  0,  their  phases  approach  the  phase  of  u»  as  |u>|  approaches 
infinity.  Thus,  in  the  limit  of  e  -  0,  the  phase  of  p,  or  R*  is 
0  for  u>  lying  on  the  real  axis  to  the  right  of  the  branch  point  lying 
on  the  positive  real  axis.  Between  the  two  branch  points  the  phase 


-8- 


is  it/ 2  and,  to  the  left  of  the  branch  point  on  the  negative  real 
axis,  the  phase  is  *.  Thus  the  phase  of  the  expression  PR'/p* 
which  appears  in  the  exponent  in  (3*2)  has  the  following  behavior 
on  the  real  axis: 

ph(PR'/u)»  0 

=  n/2 
*  0 
=  n/2 
*  It 
=  n/2 
=  n 

where 


That  the  Green’s  function  G(r-ro,u>)  may  be  considered  as 
describing  the  spatial  dependence  of  P  for  a  small  source  located 
at  rQ  may  be  readily  seen  if  one  makes  use  of  Green’s  theorem. 

One  may  show  in  general  that 

*s.->  -  4  J/  Mio'i 

-?o  4So  ’  <}-6> 


w  >  i 

(3.4a) 

A  <  w  <  Is 

(3.4b) 

A|  cos  9  |<  u)  <  A 

(3.4c) 

M  <  A  !c°?  el  . 

(3*4d) 

-A  <  a)  <  -A (cos  0  1 

( 3«4e) 

-fa  <  a)  <  -  A 

(3-4f) 

<*>  <  -  i 

(3.4g) 

Io)/|r-rJ  . 

(3-5) 

where  the  integration  is  carried  over  the  surface  of  a  small  sphere 

enclosing  the  source  (g  denoting  the  position  vector  of  points  on 

the  surface  and  denoting  the  outward-pointing  normal)  and  ^r  is 

outside  the  surface.  Then  one  may  show,  in  a  manner  similar  to  that 
( 1*>) 

outlined  by  Weston,  9  that  if  the  radius  of  the  source  is 
sufficiently  small,  the  expression  (3 becomes  approximately 


Hr,u>) 


G(r-r  ) 

71  V^O7 


IS 


7  P(r  ,w)' 
M-o  Vo 


n 


dS 


(3-7) 


I 
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The  approximation  (5*7)  is  invalid  if  either  R*  or  u  is  zero. 

In  the  limit  of  e  =  0  this  will  he  the  case  if  oj  «  A  or  w  =  A  cos  9. 

However,  we  expect  that  the  deviation  of  the  actual  pressure  field 

p(^r,t)  from  the  approximate  field  which  would  be  calculated  by  taking 

the  inverse  transform  of  (3*7)  should  be  small  if  the  source  is  small 

and  if  its  spectrum  is  not  sharply  peaked  at  these  frequencies.  With 

•  z  /2 

these  reservations,  we  shall  accept  the  product  G(r-jo,u))e  '  as 
describing  the  spatial  dependence  of  the  Fourier  transform  of  the 
pressure  field. 

The  acceptance  of  (>.7)  makes  it  possible  to  describe  the  field 
if  the  gradient  of  the  pressure  is  known  over  the  surface  of  a  small 
sphere  enclosing  the  source.  If  one  were  instead  to  consider  his 
input  data  as  consisting  of  the  recording  p^(t)  of  the  pressure  vs. 
time  at  some  point  close  to  the  source,  the  quantity  P(^r,to)  would 
be  given  by 


P(r,u>) 


G(r-r  ,w) 
G(r..-r  >w) 


(3.8) 


where  ^(w)  is  the  transform  of  p^t). 

The  Green*  s  function  G(<r-jfo,w)  has  some  interesting  features 
which  should  be  noted.  As  indicated  in  equations  (3**0,  it  corresponds 
to  a  propagating  wave  only  if  |u>|  >  ^  or  A|cos  9|<|u>|<  A.  Thus  there 
are  two  distinct  pass  bands,  the  width  of  the  lower  band  depending 
on  the  angle  0.  If  one  is  directly  above  or  below  the  source,  0  will 
be  zero  or  n  and  the  lower  band  will  have  zero  width — in  accordance 
with  the  results  derived  by  Lamb^1^  for  the  vertical  propagation  of 
plane  waves  in  the  atmosphere. 

One  may  derive  a  phase  velocity  and  a  group  velocity  for 
frequencies  lying  in  the  two  pass  bands.  These  velocities  (in  units 
of  c  and  letting  e  «  0)  are  given  by 


v 

P 


<*W(pO 


(3-9) 
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and 


[iLgAi]"1 


—  JSHlIsl.  .. 

[u  +A  (i-A2)sin  9] 


(3-10) 


where 


C  “  R'/lr-J^I  *  [n2  -  A2co820]^  .  (3-11) 


Thus  the  phase  and  group  velocities  will  depend  on  the  angle  0  in 
addition  to  the  frequency  u>.  The  phase  velocity  represents  the 
apparent  direction  with  which  wave  crests  move  in  the  radial 
direction  from  the  source,  while  the  group  velocity  corresponds  to 
the  average  velocity  with  which  a  signal  of  frequency  w  will  have 
apparently  traveled  from  the  source  to  an  observer  along  a  radial 
line  from  the  source.  One  may  readily  show  from  (3*9)  that  the 
phase  velocity  decreases  monotonically  from  •  to  0  as  increases 
from  A |  cos  9|  to  A  and  that  it  decreases  monotonically  from  ®  to  1 
as  it)  increases  from  §  to  ®.  Similarly,  one  may  show  that  the  group 
velocity  is  zero  at  w  ■  A|cos  0|  and  at  u>  a  A  and  that  it  is  positive 
and  less  than  1  for  intermediate  frequencies.  The  group  velocity 
increases  monotonically  from  0  to  1  as  the  frequency  increases  from 

b  to  00 •  In  Fig.  1  we  plot  v  and  v  vs.  for  0  -  60°. 

P  g 

In  the  frequency  bands  for  which  0  <|w|  <A|cos  0|  or  A  <|w|  < 
there  is  no  propagation  and  the  Green1 s  function  decreases  exponentially 
with  increasing  radial  distance.  The  coefficient  of  attenuation  may 
be  taken  as 


a  -  -i3C/n  •  (3-12) 

The  parameter  a  is  (^)|cos  6 1  at  u>  »  0  and  decreases  monotonically 

to  zero  as  w  increases  from  0  to  A|cos  0|.  As  u)  increases  from  A  to 
a  decreases  monotonically  from  »  to  0.  In  Fig.  2  we  plot  a  vs.  u> 
for  9  ■  60°. 
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The  factor  of  l/R*  in  the  Green*  s  function  (3*2)  will  cause  the 
Green*  s  function  to  be  singular  in  the  limit  of  e  *  0  if  |u>|  <  A  and 
|cos  6|  *  | u) | /A .  It  is  at  such  angles,  however,  that  the  approximation 
(3*7)  ceases  to  be  valid.  Nevertheless,  we  may  expect  that  the 
Fourier  transform  of  the  pressure  field  computed  via  (3*6)  for  a 
small  source  will  be  large  at  such  frequencies,  if  not  singular. 
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r i.  FORMAL  SOLOTIOW  WHEN  SOURCE  IS  ABOVE  A  FLAT  GROUND 

Ltt  us  assume  that  the  source  is  located  at  a  distance  h  (in  units 
of  H)  shore  the  ground.  The  coordinate  system  Is  chosen  such  that  the 
origin  lies  on  the  ground  directly  belov  the  source. 

To  find  the  spatial  dependence  of  the  field  corresponding  to  a 
given  frequency,  one  looks  for  a  Green*  s  function  Oh(r,w)  which 
satisfies  (5*1)  vithjc^  *  e^h,  which  satisfies  the  boundary  condition 
(2.15)  at  z  ■  0,  and  which  satisfies  the  causality  requirement  for  all 
r  for  which  z  >  0. 

The  method  of  finding  such  a  Green's  function  is  well-known. 

One  expresses  G^  as  a  Fourier-Bessel  transform, 

00 

Oj 1Ck>u)  *  J  HQ(kr)  z(k,w,z)  kdk  ,  (4.1) 

where  HQ(kr)  is  the  Hankel  function  of  the  first  kind  (defined  such 
that  its  branch  line  extends  vertically  downwards  in  the  complex 
k -plane),  and  r  represents  the  radial  distance  in  cylindrical  coor¬ 
dinates.  A  differential  equation  for  the  quantity  Z(k,u,z)  together 
with  sufficient  boundary  conditions  may  be  obtained  from  the  required 
properties  of  G^.  In  this  manner,  one  may  obtain  the  following 
expression  for  Z(k,u>,z): 

.  K  [.‘♦Ml  . -pL.it l»»li ,  (k.a) 

2CT  L  B  +it  J 


♦  -  [p2  -  iA2/fl2]*  (4.3) 


where 


The  phase  of  t  Is  chosen  such  that;  If  k  Is  real,  t  has  no  branch 
line  In  the  upper  half  of  the  complex  w-plane  and  such  that  the  phase 
of  t  approaches  the  phase  of  w  as  |u)|  approaches  infinity  if  Im(ci>)  >  0. 
(This  will  insure  that  conforms  to  causality  in  the  sense  previously 
described.) 

The  phase  requirements  on  i  in  the  u> -plane  enable  one  to  derive 
the  phase  of  f  for  all  real  k  for  any  given  real  u.  In  general,  for 
<  not  identically  zero,  one  may  specify  the  phase  of  t  uniquely  for 
all  real  k  and  real  u  by  requiring  the  phase  of  f  to  be  between  0  and 
x.  The  specification  may  be  extended  to  include  complex  values  of  k 
by  placing  branch  lines  in  the  k -plane  and  requiring  that  the  phase 
of  t  be  continuous  everywhere  except  at  these  branch  lines.  For  «  >  0, 
the  placing  of  these  branch  lines  is  restricted  by  the  fact  that  no 
branch  line  may  cross  the  real  axis,  since,  in  this  case,  the  phase 
must  be  continuous  along  the  real  axis.  The  simplest  placing  is  to 
take  the  branch  lines  as  extending  vertically  upwards  from  branch 
points  lying  above  the  real  axis  and  vertically  downwards  from  those 
lying  below  the  real  axis.  One  may  readily  show  that  the  branch 
points  of  f  in  the  k-plane  lie  in  the  first  and  third  quadrants  if 
(4  >  0  and  in  the  second  and  fourth  quadrants  if  u>  <  0.  Thus,  in  the 
limit  of  «  ■  0,  branch  lines  extend  upwards  from  branch  points  lying 
on  the  positive  real  axis  and  downwards  from  branch  points  lying  on 
the  negative  real  axis  if  u>  >  0,  while  the  converse  is  true  if  u>  <  0. 

The  specification  of  the  phase  of  ♦  as  outlined  above  reduces 
to  the  following  for  c  *  0  and  real  k: 


1. 

For  |W|  >£  , 

Ph(t)  -  Ph(u) 

w 

1  <  |fWn| 

(4.1*) 

-  */2 

w 

>  lpU)Al| 

(4.4b) 

2. 

For  £  >  |w|  >  A, 

Ph(0  -  x/2 

all  k 

(*•5) 

3- 

For  |  ut  |  <  A, 

Ph(t)  -  x/2 

w 

<  IfVnl 

(4.6a) 

-  x  +  Ph(u) 

w 

>  |Pw/n| 

(4.6b) 
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The  relation  given  by  (4,6b)  should  particularly  be  noted.  More 
naive  considerations  might  lead  one  to  believe  that  the  phase  of  '!/ 
is  either  equal  to  the  phase  of  w  or  to  jt/2  for  real  k. 

For  e  >  0,  the  integration  in  (4.1)  may  be  taken  as  proceeding 
along  the  real  axis.  The  poles  of  Z(k,u>,z)  will  lie  off  the  axis-- 
specifically,  at  the  points 

k  =  ±  (u>  +  U/2)  (4.7) 

for  sufficiently  small  e .  If  e  =  0,  the  poles  will  lie  on  the  real 
axis,  and  the  contour  must  pass  below  the  pole  on  the  positive  real 
axis  and  above  that  on  the  negative  real  axis  if  u>  >  0.  If  w  <  0, 
the  converse  is  true. 

Let  us  now  note  that  the  free-space  Green's  function  G^-^hjw) 
given  by  (3*2)  may  also  be  expanded  in  the  form  (4.1).  The  appropriate 
expression  for  Z(k,u),z)  may  be  readily  shown  to  be 

-itt-  ei’Hz-hl  (4.8) 

2n  t 

in  a  manner  similar  to  that  outlined  above.  Thus,  with  the  decompo¬ 
sition  of  Z  indicated  by  (4.2),  our  Green's  function  G^  appears  as  a 
sura  of  three  terms: 

Gh(£,u)  =  Gtazh'u)  +  °(Sf'iSzh*u)  +  I(r,z+h,w).  (4.9) 


The  first  term  represents  the  free-space  Green' s  function  for  a 

source  located  at  r  »  e  h,  while  the  second  term  represents  the 
**o  **z 

free-space  Green's  function  for  the  image  source.  The  third  term 
is  given  by  the  integral 


l(r,z+h,w) 


«  2 

f  —I — 


if(z+h) 

e  HQ(kr)  kdk  . 


(4.10) 
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This  term  is  present  as  a  direct  result  of  the  fact  that  dP/dz  is  not 
identically  zero  at  the  surface  of  the  earth. 

We  may  regard  the  product 

*Z/2  Gh(£'u) 

as  describing  the  spatial  dependence  of  the  pressure  field  from  a 
harmonic  source  or  as  describing  the  spatial  dependence  of  the 
Fourier  transform  of  the  field.  This  follows  from  Eqn.  (3.8). 


V*  THE  FIELD  BELOW  THE  SOURCE 


The  fictional  parameter  c  is  of  no  further  use  to  ue  and  will 
henceforth  he  set  to  zero.  Also,  since 

°h(r,  -«)  -  -  0 *(r,w)  (5.1) 

for  real  w,  we  shall  limit  ourselves  to  positive  w. 

Since  f  is  even  in  k,  the  Integral  in  (4.10)  may  also  be  taken 
in  the  fora 

-iu  r°  2B2  *♦(*•*) 

l(r,z+h,u>)  -  — §  |  — =5 -  e  J  (kr)  kdk  (5-2) 

u  Jo  t(B +!♦)  ° 

with  the  contour  passing  below  the  pole  at  k  *  w.  This  fora  is  most 
convenient  for  studying  the  field  at  small  r.  If  we  take  r  ■  0,  the 
Bessel  function  becomes  1  and  the  variable  of  Integration  may  be 
changed  to 


s  -  -  (B2+it)(z+h) 


to  give 

2  <*> 

l(0,z+h,oj)  -  (2B2A0  e"B  f  (e~a/s)  ds  ,  (5*5) 

J-s 

o 

with 

•0  -  (B2  +  10)(z+h)  ,  (5^) 

and  the  contour  passing  below  the  pole  at  s  •  0.  The  integral  in 
(5.5)  is  the  exponential  Integral  of  complex  argument  and  has  been 
tabulated  by  the  National  Bureau  of  Standards If  | »Q |  is  very 
large,  the  Integral  is  very  nearly  -exp(sQ)/so  and  we  have 
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l(0,z+h,u) 


-2B 


>i(B2+ip) 


1P( s+h) 
(s+h) 


(5-5) 


If,  on  the  other  hand,  |*o|  is  email. 


l(0,z+h,w) 


-B2(z+b) 
e  loge 


{-(B2+ip)(  s+h)e*5772} 


(5-6) 


with  the  branch  cut  for  the  logarithm  on  the  positive  Imaginary  axis* 
We  note  that  when  (5**0  Is  applicable  the  ratio  of  I  to 
Gfrff^hjO))  at  r  «  0  is 

-2B2  (5.7) 

B2  +  ip 


This  ratio  Is  vanishingly  small  for  large  frequencies!  but  may  be 
very  large  If  w  is  close  to  A*  If  (5*6)  is  applicable,  the  ratio 
Is 


2B2(z+h)  loge{-(B2+ip)(s+h)e5772}  ,  (5.8) 

which  le  singular  for  w  ■  A. 

Let  us  now  consider  (5*2)  when  r  Is  not  Identically  zero.  If 
u  <  A  and 

r  <  (z+h)tan  $c  ,  (5*9) 


where 


tan  9 

c 


[<*A.>2  -  l]*-. 


(5.10) 
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the  contour  in  (5*2)  may  be  deformed  to  the  negative  imaginary 
With  an  appropriate  change  of  variable,  (5*2)  then  becomes 

-ilu.1  i“  -|  +  |(z+h) 

l(r,z+h,w)  =  I  f(%)  I  (?r)  e  5d5 

0)  J  o  c 

where  IQ(?r)  is  the  Bessel  function  of  imaginary  argument, 

1*1  ■  [lel 2  *  (WA-)8]4  . 

and 

t(%)  =  — 

|*|(|t|-B ) 

is  a  positive  real  function  of  * .  It  is  immediately  clear  that  the 

magnitude  of  (5*11)  increases  with  increasing  r  and  decreasing  z+h. 

It  is  also  clear  that  the  phase  of  I  will  not  depend  on  r  or  z. 

Let  us  note  that  the  latter  property  is  shared  by  G^r+e^h) 

if  a)  <  A  and  (5*9)  is  satisfied,  and  that  it  also  holds  for  GQr-e^h) 

if  r  <  lz-h|tan  9  .  We  may  accordingly  conclude  that 
*  c 


axis. 

(5-n) 


s  -  -  “hi0!.’  <5a2> 

is  zero  everywhere  above  and  below  the  source  within  a  cone  of  apex 
angle  9c  with  the  apex  at  the  source. 

The  quantity  introduced  above  gives  the  spatial  dependence  of 
the  acoustic  intensity  for  a  harmonic  source.  (This  may  be  seen  by 
taking  the  product  of  the  pressure,  as  given  by  (2.l6),  and  the 
complex  conjugate  of  the  particle  velocity,  whose  transform  is  given 
by  (2.13),  and  then  integrating  over  time.)  Thus  we  may  conclude  that 
for  a)  <  A,  there  is  no  energy  propagation  in  the  cone  of  apex  angle  0c 
Let  us  now  consider  (5*2)  in  the  case  where  u>  is  between  A  and^  . 
In  this  case  we  may  consider  the  integral  as  being  a  sum  of  its 
principal  value  plus  an  integral  around  a  small  semicircle  passing 
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belov  the  pole.  The  latter  may  be  evaluated  by  the  method  of  residues. 
The  contribution  to  I  from  the  principal  value  of  the  integral  will  be 
a  real  number  as  t  has  a  phase  of  n/2  along  the  real  axis.  Thus  the 
imaginary  part  of  I  will  be  given  by  the  contribution  from  the  pole: 


lm(l) 


-B2(z+h) 

e  Jo(u>r). 


(5.13) 


This  also  gives  the  Imaginary  part  of  G.  (r,u>)  . 

n 

The  quantity  l(r,z+h,u>)  is  related  to  the  Green's  function 
G(£tieQi,u)  by  the  equation 

(d/dz  +  B2)  I  -  -  2B20  (5. 1*0 


as  may  be  seen  from  an  examination  of  (4.10).  This  equation  is  useful 
if  one  seeks  to  derive  an  expression  for  the  acoustic  intensity  jS^. 

The  z-coraponent  of  E^may  be  found  directly  from  (5. 12)  if  one  uses 
the  above  relation  and  the  fact  tnat  the  two  Green's  functions  are 
real.  In  this  manner,  one  obtains 

S2  «  -(2/n2)  lm(l)  {(|j  -  B2)G(r+eh,u>)  +  (fj  *  B2)G(r+<>zh,U)}  . 

(5.15) 

The  r-component  may  be  found  if  one  uses  the  fact  that  the  divergence 
of  S  is  zero.  Thus 


-<l/r)  f 


(3S  /dz)  rdr. 
z 


(5.16) 


After  some  algebraic  manipulations  and  an  integration  by  parts,  this 
becomes 
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Sr  =  -(2/u,2)  {Ml)  Ij  [G(^zh) +  G(£-Szh] 

-  [G(,^-ezh)  +  °Cr-«zh)]  §7  Im(1)} 

+(4B2n‘2/r)  j  Im( I)  (B2  -  |j)G(r+Szh)  rdr  •  (5*17) 


From  (5*15)  and  (5*17)  it  is  clear  that  the  intensity  near  the 
source  is  of  the  form 


S 


,  2  (m  h)  e  *  n 

4nB  **  z 

“T£~  2  2  ^  77^13/2 

*  '.or  +  ^  (z-h)  " 


(5-18) 


which  represents  a  flow  of  energy  radially  outward  from  the  source* 

2  2 

However,  since  ^  ^  w  ,  the  flow  is  not  spherically  symmetric. 

More  energy  flows  out  vertically  than  horizontally.  We  also  note 
that  as  much  energy  flows  out  in  the  upward  direction  as  flows  out 
in  the  downward  direction. 

One  interesting  feature  of  our  expressions  for  S  is  that  S 

**•«*>  z 

is  proportional  to  the  Bessel  function  JQ(u>r).  Whenever  r  is  such 
that  wr  is  a  root  of  the  Bessel  function,  S  will  be  zero  and  the 
energy  flow  will  be  horizontal  for  all  points  on  the  cylinder  of 
radius  r. 

In  the  limit  of  large  r  the  intensity  is  entirely  in  the 
horizontal  direction  and  is  given  by  the  last  term  of  (5-17): 

2  oo 

Sr  =  (8kbV^)  e"B  (z+h^r~1f  JQ(wr)(B2  -  b/bz.)G(r+^h,u)  rdr  . 

(5-19) 

This  indicates  that  the  intensity  decreases  inversely  with  r  at 
large  distances. 


-21- 
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The  nature  of  the  transition  from  (5*18)  to  (5*19)  is  best 
illustrated  by  means  of  a  numerical  example.  We  choose  u>  *  .W3 
and  h  *  5*0*  (in  conventional  units  h  would  be  approximately  ^0  km.) 
The  results  of  the  calculation  are  given  in  Pig.  3.  There  we  show 
the  typical  paths  along  which  energy  would  flow  away  from  the  source. 
(These  lines  may  be  considered  as  being  the  sound  rays.)  The  figure 
also  shows  the  variation  of  the  magnitude  of  the  intensity  with 
distance.  The  downward  bending  of  the  rays  in  the  region  above  the 
source  is  reminiscent  of  a  water  fountain  and  reminds  us  that  we  are 
considering  a  problem  in  which  the  earth’s  gravity  is  important. 
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VI  .  APPLICATION  OF  THE  METHOD  OF  STEEPEST  DESCENTS 


The  case  when  wr  is  large  is  best  studied  within  the  context 
of  the  method  of  steepest  descents.  If  one  follows  the  rationale  of 
this  method,  the  Hankel  function  in  (4.10)  is  replaced  by  its  asymp¬ 
totic  expression  and  the  contour  is  deformed  to  one  which  passes 
through  the  saddle  point  of 


eikr+it(z+h)  (6.1) 

and  which  goes  along  the  path  for  which  this  quantity  decreases  most 
rapidly  with  increasing  distance  from  the  saddle  point. 

The  saddle  point  is  located  at 

k  «  k 

s 


where 

tan  9  ■  r/(fc+h)  ,  (6*3) 

Q  M  (u)2  -  A2cos20)^  (6.4) 

and  9  is  between  0  and  «/2.  The  path  of  steepest  descents  is  given  by 

,  (k  -k_)(k  -|k  |) 

kT  “  (wA0  r  "  ’5  ■’  tan  0 >  u  £  (6.5a) 

[<Vk.>  *«] 

-  (w/ft)(k2  ♦  Q2)^,  #  >  w  >  A  (6*5b) 

*  (w/M)  R  -£ — -^r-  tan  ft,  A  >  w  >  A  co.  ft  (6*5c) 

oh,  -VK  -V4 


(6.2) 
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and 


-  0  , 

0  <  u)  <  A  cos  0  , 

(6.5d) 

= 

lK/(u  Sin  0)1  , 

(6.6a) 

Q  = 

!(u>P  cos  9)/c|  , 

(6.6b) 

kN  ' 

lksl  -  Q  , 

(6.6c) 

^  “ 

|kg|  +  <a , 

(6.6d) 

and  k  and  k  denote  the  real  and  imaginary  parts  of  k,  and,  in  (6.5c), 
K  X 

kR  is  understood  to  be  between  k^  and  k^.  If  we  let 

kj  » 

represent  the  distance  of  the  branch  points  from  the  origin,  then  it 
is  readily  verified  that  kg>  k^>  Jkg|  for  w  >  £,  |kg|  >  k^  for 
£  >  w  >  A,  and  k^  >  |kg|  >^>^>^>0  for  A  >  w  >  A  cos  0. 

In  Fig.  4  we  show  sketches  of  the  path  of  steepest  descents  for  the 
various  ranges  of  w. 

The  deformation  of  the  contour  in  (4.10)  to  the  path  of  steepest 

descents  is  valid  for  w  >  £  only  if  k^>  uj,  and  is  valid  for 

A  cos  0  <  u)  <  A  only  if  kg<  w.  If  this  is  not  true,  one  must  add 

the  residue  from  the  pole  at  k  «  u>  to  the  integral  I  along  the 

SD 

path  of  steepest  descents  to  obtain  (4.10).  This  residue  is 


-24- 


Res 


2*i  -B  (z+h) 


H  (u)r) 
ov  ' 


2 

\ro  //  -B  (z+h)  iwr 

«  e  '  (2B  /^)[2n/((or)r  e  v  e 


(6.7a) 


(6.7b) 


The  condition  k^,  >  u>  for  u>  >  £  is  equivalent  to  to  >  u>^(  0), 
while  the  condition  k^,<  u)  for  A  cos  8  <  w  <  A  is  equivalent  to 
u>  <  (0),  where 

(2  cos20)^  2(0)  *  \  +  A2cos20  +  j(£  -  A2)2  +  A2(i  -  A2)sin220p  . 

(6.8) 

(It  is  readily  seen  that  increases  monotonically  f rom  ^  to  00 
while  decreases  monotonically  from  A  to  0  as  0  goes  from  0  to  n/2.) 
The  integral  (4.10)  is  thus 


I(r,z  +  h,u>)  =  ISD  +  U(w,0)  IReg  (6.9) 

where  u(w,0)  is  1  for  u>^(0)  >  w  >  uv>(0)  and  is  otherwise  zero. 

The  integral  along  the  path  of  steepest  descents  may  be  evaluated 
approximately  in  the  usual  manner,  giving 

ISj)  «  [t(w,9)  -  l]  G(r+e2h,co)  ,  (6.10) 

where 

T(o>,e)  *  sss. 2-2?  .  (6.11) 

i(|Ap/C)  cos  0  +B 


It  is  clear  that  |t|  is  1  if  w  >  J  or  A  cos  0  <  w  <  A  and  that  t  is 
real  and  less  than  -1  if  £  >  u>  >  A.  Also,  t  is  real  and  greater  than 
1  if  a)  <  A  cos  0.  In  Fig.  5,  we  plot  the  phase  and  magnitude  of  t 
versus  a)  for  several  representative  values  of  0. 
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With  the  approximation  (6.10),  the  field  G^frjZ,^)  becomes 

G.  ( r,z,io)  =  G(r-e  h,io)  +  TG(r+e  h,w)  +  U  1^,  •  (6.12) 

h  \v  Res 

The  three  terms  in  this  equation  may  be  conveniently  labelled  as  a 
direct  wave,  a  reflected  wave,  and  a  surface  wave,  respectively. 

The  parameter  T  thus  corresponds  to  a  reflection  coefficient. 

Let  us  now  study  the  behavior  of  (6.12)  for  given  u>  and  h  vs. 
r  for  z  =  0.  In  this  event  the  two  Green's  functions  are  equal. 

Also,  if  w  >  J  or  u  <  A,  the  presence  of  the  factor  U  in  (6.12) 
implies  that  the  surface  wave  will  nob  be  present  unless 

r  >  rQ  =  |3ti/(B2(o)|h  .  (6.15) 

It  need  scarcely  be  pointed  out  that  this  parameter  rQ  has  no 
physical  significance.  It  is  unlikely  that  the  field  should  change 
abruptly  at  r  =  rQ.  Furthermore,  at  this  value  of  r,  the  surface 
wave  may  be  of  insignificant  value  compared  to  the  sum  of  the  direct 
and  reflected  waves  if  the  conditions  for  the  validity  of  the 
approximation  (6.10)  are  satisfied.  Nevertheless,  the  theory  does 
give  us  a  natural  small  distance  cut-off  for  the  surface  wave. 

This  indicates  that  Eqn.  (6.12)  will  be  qualitatively  correct  for 
all  r  for  frequencies  in  these  two  bands  if  the  height  of  the  source 
is  sufficiently  large. 

As  r  increases  from  r  ,  the  relative  contribution  of  the  surface 

o 

wave  to  G,  will  increase.  This  follows  for  two  reasons.  First,  the 
h 

surface  wave  falls  off  more  slowly  with  r  than  the  direct  and  reflected 
waves.  Secondly,  the  direct  and  reflected  waves  will  tend  to  cancel 
each  other  at  large  r. 

A  characteristic  distance  r^  may  be  defined  as  that  value  of  r 
for  which 


l^esl  “  I1  +  T>  !Gtezh'u,)!  • 
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For  values  of  r  greater  than  the  surface  wave  will  be  dominant, 

while  the  converse  will  be  true  if  r  <*  r,  .  For  go  >  h  or  to  <  A  this 
parameter  may  be  shown  to  be  approximately 


r 


1 


h2 


(6.1k) 


2  2 

provided  this  expression  is  much  greater  than  |u|  h/(B  u>) 
or  ju|h/co.  (These  criteria  will  be  met  if  h  is  large.)  In  the  high 

frequency  limit,  r  is  proportional  to  co,  while,  in  the  low  frequency 

1  -1 
limit,  it  is  proportional  to  <i>  ♦ 

If  A  <  cj  <  Eqn.  (6.12)  gives  no  small  distance  cut-off  to  the 

surface  wave.  Although  the  expression  is  not  valid  for  r  identically 

zero,  we  expect  it  to  be  approximately  valid  for  any  r  for  which 

tor  >  1,  assuming  h  is  large.  This  is  confirmed  by  Fig.  4,  which 

shows  a  surface  wave  present  for  almost  all  points  on  the  ground 

except  directly  below  the  source.  For  r  =  0,  the  field  may  be  found 

by  use  of  (5*3)*  In  the  limit  of  large  h,  (5*5)  is  applicable  and 

we  find 


G(  0,0, u>) 


-g|3l  e~lglh 
B2  -  |0| 


(6.15) 


It  is  interesting  that  this  has  sign  opposite  to  that  of  the  free 
space  Green* s  function. 
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VII.  DISCUSSION  AND  CONCLUSIONS 


The  explicit  solution,  Eqn.  (3*2),  for  the  case  when  the 
presence  of  the  ground  is  ignored,  represents  the  correct  solution 
of  the  problem  originally  considered  by  Sretenskii The  method 
of  deriving  this  solution  (i.e.,  of  proceeding  from  the  equations  of 
hydrodynamics  in  a  linearized  form,  then  seeking  a  Green’s  function 
conforming  to  causality,  and  finally  using  Green's  theorem  to  find 
the  field  in  the  limit  of  an  infinitesimal  source)  appears  to  be 
the  most  consistent  approach,  and  it  is  gratifying  that  our  solution 
has  such  a  relatively  simple  form.  The  prediction  of  the  two  non¬ 
propagating  frequency  bands  is  perhaps  the  most  interesting  consequence 
of  our  theory.  It  appears  difficult  to  give  a  simple  qualitative 

explanation  for  the  presence  of  these  bands.  We  may  mention,  however, 

(17) 

that  the  frequency  w  *  A  is  Brunt's'  resonant  frequency  for  the 
isothermal  atmosphere.  The  frequency  w  ■  J  is  also  a  characteristic 
frequency  for  the  isothermal  atmosphere  and  appears  in  Lamb's^1^ 
theory  as  the  dividing  point  between  propagation  and  attenuation  of 
vertically  propagating  plane  waves. 

When  the  presence  of  the  ground  was  incorporated  into  the  theory, 
the  solution  to  the  nroblem  became  more  complicated.  Although  it 
could  be  given  in  integral  form,  the  term  l(r,z+h,u>)  could  not  be 
evaluated  exactly.  We  were  able,  however,  to  derive  an  explicit 
expression  for  the  case  r  ■  0,  enabling  one  to  compute  the  field 
anywhere  on  the  vertical  line  passing  through  the  source.  We  were 
also  able  to  give  the  solution  in  the  limit  when  the  method  of 
steepest  descents  was  applicable.  This  led  to  a  representation  of 
the  solution  as  a  sum  of  a  direct  wave,  a  reflected  wave  (with 
reflection  coefficient  (6,11)),  and  a  surface  wave.  The  surface  wave 
is,  of  course,  the  counterpart  of  the  gravity  modes  discovered  in 
other  models  of  the  atmosphere. 

For  an  understanding  of  the  nature  of  the  field  when  w  <  the 
study  of  the  acoustic  intensity^  was  useful.  We  found  that ^ was 
zero  if  a)  <  A  everywhere  above  and  below  the  source  within  a  cone  of 
apex  angle  0c*  This  was  just  what  might  have  been  expected  after 


-28. 


considering  the  free-space  problem. 

The  field  for  frequencies  in  the  band  A  <  u>  <  £  has  some 
interesting  properties  which  were  brought  out  by  our  theory. 

In  the  free-space  problem,  no  propagation  would  take  place  and 
S  would  be  everywhere  zero.  However,  the  presence  of  the  ground 
permitted  the  existence  of  a  propagating  surface  wave.  We  were 
able  to  derive  relatively  simple  expressions  for  the  intensity 
above  the  ground.  A  numerical  example  (Fig.  3)  showed  that  the 
energy  flowing  out  of  the  source  had  an  overwhelming  tendency  to 
concentrate  in  the  lower  layers  of  the  atmosphere. 

For  arbitrary  frequencies  and  at  sufficiently  large  distances, 
the  field  near  the  ground  would  be  entirely  given  by  the  surface 
wave.  However,  for  high  frequencies  or  very  low  frequencies,  the 
surface  wave  would  not  predominate  unless  the  distance  were  very 
large . 
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A  cos  9  A  0.5 

Angular  frequency  in  units  of  C/H 


Fig.  1.  —  Phase  velocity  V_  and  group  velocity  V„  (units  of  c) 

r  Gr 

versus  u>  (units  of  c/H)for  free-space  propagation 
from  a  point  source  along  a  line  making  an  angle  of 
6o°  with  the  vertical. 


A  0.5 


A  cos  9 

Angular  frequency  u>  in  units  of  C/H 

Fig*  2.  --  Attenuation  coefficient  a  (unit 6  of  H  versus 
w  (units  of  c/h)  for  free-space  propagation  from 
a  point  source  along  a  line  making  an  angle  of 
6o°  with  the  vertical. 


-31- 


H  jo  sjjun  uj  iq6i9H 


Fig.  4.  —  Sketches  of  the  cooplex  K -plane  showing  relative  location 
of  branch  lines,  poles,  and  path  of  steepest  descents  for 
various  ranges  of  the  angular  frequency  u.  (The  saddle 
point  Is  denoted  by  Kg«) 


Fig.  4.  —  Sketches  of  the  complex  K-plane  shoving  relative  location 
of  branch  lines,  poles,  and  path  of  steepest  descents  for 
various  ranges  of  the  angular  frequency  u>.  (The  saddle 
point  is  denoted  by  Kg*) 
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